Exact Langevin equation for electrons in thermostat, 
Fokker-Planck kinetics and flicker noise 



Yuriy E. Kuzovlev 
Donetsk Physics- Technical Institute of NASU, Donetsk, Ukraine 

(February 2001) 



Abstract 

Generalized Langevin equation for characteristic functional of many-electron 
system dynamically interacting with a thermostat and besides subjected to 
external perturbation and observation is derived and formulated in terms of 
one-particle but stochastic density matrix. The scheme naturally includes 
also direct electron-electron interactions. Corresponding approximate Fokker- 
Plank kinetic equation is discussed. We point out that it possesses continuous 
degeneracy of stationary solution and therefore can hide long-correlated fluc- 
tuations of electron transport even if thermostat noise is short-correlated. 

A. Introduction 

The low- frequency flicker noise (1/f noise) is of growing interest among scientists. Usually, 
investigators interpret it as manifestation of some slow kinetics, for instance, of thermally 
activated switchings with a broad variety of memory-times (life-times, relaxation times, 
etc.). From this point of view, 1/f-noise says that present events in a system depend on 
its old history. Opposite principal possibility is that 1/f-noise is due to forgetting history, 
concretely, forgetting how many events took place in the past, which means that a system 
has no reasons to keep well certain "number of events per unit time". Then 1/f spectrum 
reflects only statistical freedom of flow of events but not real cause- consequence connections 
between them. From this point of view, 1/f noise in electronics and, say, 1/f fluctuations of 
flowing of the river Niel are similar phenomena. General ideology of this approach and its 
concrete realizations were published elsewhere [1-11]- 

The two different ideas imply different programs for future theory. In our approach, one 
should reconsider derivations of kinetic models of dissipation (irreversibility) and noise from 
underlying Hamiltonian dynamics, avoiding the ansatz that "numbers of events per unit 
time" (i.e. collision integrals, etc.) have quite definite (and well achievable) limits as being 
averaged over time. The suspicion that such quantities may be uncertain was claimed and 
argued by Krylov [12] although with no relation to 1/f-noise. It must be emphasized that 
violation of ergodicity with respect to events, i.e. transitions between phase space regions, 
does not mean violation of ergodicity with respect to visiting these regions, therefore, usual 
ergodic theorems are not under doubts. Also we have no doubts of importance of slow 
kinetics (thermally activated or other). 



In particular, it is interesting to elucidate whether many- electron system interacting 
with a thermostat can possess 1/f type or other flicker fluctuations in electron transport if 
autonomous thermostat dynamics itself is free of such fluctuations. Any definite answer for 
this question would be equally useful. In this work we make several steps in this direction 
and argue that the answer is perhaply positive. 

The object under our interest is characteristic functional of electron current. First, we 
derive formally exact "Langevin equation" for fermions subjected to Hamiltonian interaction 
with a thermostat, to external dynamical perturbation (bias electric field) and besides to 
external observation. Secondly, we show that operator-valued thermostat variables in this 
equation can be effectively replaced by doubled c- number- valued variables (" noises" ) . As the 
consequence, the equation allows to perform evaluation of fermionic statistical sums before 
averaging with respect to thermostat. Therefore, the problem reduces to analysing statis- 
tical properties of one-electron but instead stochastic ("twice randomized") density matrix. 
Next, supposing that electron-thermostat interaction is weak and short-correlated, we de- 
rive approximate Fokker-Planck equation for probability distribution of the density matrix, 
and then consider corresponding kinetic equation for its eigenvalues which just determine 
the characteristic functional. Finally, we attract attention to that stationary solution of 
the Fokker-Planck equation is continuously degenerated, and that this fact can signify exis- 
tence of long-correlated fluctuations in electron transport (flicker type excess current noise) . 
Mathematically rigorous proof (or refutation) of this statement needs in analysis of spectra 
of comphcated non-self- adjoint second-order differential operators. This will be next stage 
of the work. 



B. Model 

Let a quantum system consists of an essential "dynamical subsystem", with Hamiltonian 
He{t) , and a "thermal bath" ("thermostat"), with Hamiltonian Hb . We assume that (i) 
quantum states of these two parts belong to different independent Hilbert spaces, and (ii) 
their interaction is described by Hamiltonian 

H,^,^J2Bj*D, (1) 
j 

where Dj and Bj act in dynamical space and bath space, respectively. This is the case in 
most of physically meaningful models. For concreteness, we may take in mind that dynamical 
subsystem (to be abbreviated "DS") is formed by electrons obeying secondary quantization, 
while bath by bosons (for instance, phonons). Thus the total Hamiltonian is 

H{t)^ Heit)+Hb + Hint (2) 

Besides the interaction, DS may be subjected to an external quasi-classical (c-number) 
perturbation, u{t) . In many cases it also can be described by bilinear form, 

He{t) = He- u{t)Q , (3) 

with Hf. being unperturbed Hamiltonian and Q a conjugated quantum variable. For exam- 
ple, if u{t) is an external voltage then Q represents amount of charge transported by free 
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(conducting) electrons, or if u{t) is applied electric field then Q is polarization. Generally 
speaking, operator Q docs not have definite expression, but the corresponding flow (current), 
i[H{t),Q] , is well defined (square brackets will serve for operator commutation, [^, i^] — 
XY — YX ). We assume that (iii) Q belongs to the dynamical subspace (or to some its 
extension) and (iv) Q commutes with all the operators Dj , and thus [Hint, Q] — . This 
means that thermal bath prodicing dissipation and decoherence in DS at the same time do 
not contribute to the current in a direct way, as usually in models of electron transport. 
Then the current operator 

J^i[H{t),Q] ^i[H„Q] (4) 

also purely belongs to the DS subspace, and 

[H,, J] = (5) 

Of course, in any non-trivial case operators Dj do not commute with free DS Hamiltonian, 
[He, Dj] ^ , and corespondingly [J, ^ , although commutation [i^e, J] = is possible. 



C. Quantum characteristic functional 



The things under our interest will be fluctuations in DS, firstly in the current, J , and 
its response to external perturbation. To find these, one should consider time-dependent 
operators in the Heisenberg representation. 



Jit) = S+(t)JS(t) , S(t) = gxp \-i f H{t')dt' 

L Jo 



(6) 



where superscript "+" denotes Hermitian conjugation, and symbols ^xp and ex|5 serve for 
chronological and anti-chronological time-ordering exponents, respectively. In the Heisen- 
berg picture, fiuctuations are described by statistical moments (J(ti)... J(t„)) , where angle 
brackets is averaging via some statistical ensemble of initial states. 



(...) = Tr , 



(7) 



with Pin being an initial density matrix. 

To have real-valued moments, it is necessary to choose a certain symmetrization rule, in 
other words, some quantum generalization of classical characteristic functional 



exp 



v{t')J{t')dt' 



where v{t) is arbitrary test function. The first and second moments are indifferent to the 
choice. One of possible definitions is 



H(v) = 



2io 



v{t')J(t')dt' 



6xp 



t vit')Jit')dt' 
Jo 



1 rt 
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(8) 



It corresponds to following time-ordered symmetric product: 
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J{h)...J{tn) = {..{J{tl)oJ{t2)}..oJ{tn-l)}oJ{tn) if > ^2 > > , 



where the circle means Jordan symmetric product, X oY = ^{XY + YX) . 

The convenience of this variant of symmetrization is that it allows to reduce calculation 
of the integral quantity S(f ) to solution of a differential equation. With the help of well 
known general recipes for disentangling operator exponents we can rewrite expression (8) as 



E{v) 



(ex 



v(t')J + iH{t') 



dt) ^ (/^ 



v(t')J -iH{t') 



dt' 



(9) 



Because of general properties of the Trace operation this means that 

E{v) = Tr p{t) , 
where operator p{t) obeys the equation 

v{t)J o p{t) + z[p{t), H{t)] = v{t)J o p{t) + L{t)p{t) , 



d , . 



(10) 



with initial condition p(0) = pm ■ Equation (10) remains valid also in classical limit, if replace 
L(t) by Liouville operator and J by phase function of canonical variables, respectively. At 
f (t) = it turns into mere von Neumann (Liouville) evolution equation for density matrix 
(probability distribution). Hence, the testing is equivalent to non-Hamiltonian perturbation 
which can not be introduced by means of Poisson brackets. 



D. Free bath representation 



We will suppose that the only significant role of thermal bath is destroying coherence of 
DS evolution, i.e. production of noise, dephasing and dissipation. Otherwise (for instance, 
if the bath changes ground state of electronic subsystem), we must redefine division of the 
whole system into two parts. Then, it is reasonable to speak the free bath language in which 
bath behaves independently on DS, while the latter is disturbed by bath and governed by 
the Hamiltonian 



//(S,i)=i/e(i)+E^iW^.- 



(11) 



Here the operators 



Bj{t) = exp{iHut)BjCxp{-iHbt) 



describe free bath evolution. Taking into account the equality (5), it is easy to verify that 
expressions (6) and (9) are indeed identical to 



J(t) = SMJSBit) , SB(t) = ^ \-i t H(B,t')dt' 

L Jo 



(12) 



=(.) = (e-^(/;[ 



-v{t')J + iH{B,t') 



dt')^,{i; 



-v{t')J-iH{B,t') 



dt' 



(13) 



4 



E. Quantum thermostat to classical noise conversion 



The "free bath representation" (12)-(13) allows to replace quantum bath by an effective 
classical disturbation (noise), or to be precise, by pseudo-classical one. In perfect sense this 
can be made under assumption that (v) initial statistical ensemble has the factored form 



Pin ^ Pe* Pb : 



(14) 



with no initial correlations between two subsystems, and therefore the complete averaging 
exactly divides into two stages: 



(•••) = ((•••)e). , (...)e = TrJ...}p, , (...), ^ Tr,{...}p,, 

where Trg and Tih are traces in dynamical and bath subspaces. 
Consider the construction 

E{v) = ^cxi^ [it;(tO J + iHiv, O] dt'^ Sep (^J^ [lv{t')J - iH{C, t' 



(15) 



dt' 



(16) 



with different left-hand and right-hand Hamiltonians defined by replacing B{t) in (11) by 
r){t) and C,{t) , and i]j{t) and ^j{t) being formally treated as classical (c-number valued) 
noises. Obviously, this expression coinsides with (13) if interpret (...)^ as averaging with 
respect to noises r]{t) and ^{t) under an appropriate choice for their statistics. To get the 
exact identity, we should introduce their characteristic functional by 



exp J T.MtMt) + bj{t)rjj{t)]dt) 
^ i ' b 



(17) 



= Tib ex|5 



J2b,{t)B,{t)dt 



5xp 



J2aj{t)B^{t)dt 



Pb , 



where aj{t) and bj{t) are test functions. As consequence, instead of (10), wc can write 

E{v) ^{Tr, p{t)), , 

where now "density matrix" p{t) obeys the generalized von Neumann equation 
d 



dt 



p{t) = vit)Jo p{t) + L{t)p{t) , L{t)p = iipHi-n, t) - H{^, t)p} , 



(18) 



with initial condition p(0) = Pe ■ It concerns DS only and evolves purely in the dynamical 
subspacc, while ri{t) and ^(t) play a role of noise perturbation externally applied to DS. 

The ansatz (14) is natural approximation if the interaction Hint is sufficiently weak. 
But in general case we also may start from the ensemble (14), if relating E{v) to later 
times than characteristic time of mutual relaxation of subsystems, and thus considering 
their thermodynamical interaction as result of dynamical time evolution. 
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F. Generalized Langevin equation 

The convenience of the pseudo-classic representation introduced by Eqs. 16-18 is that one 
can deal with ri{t) and ^{t) as objects which commute with both themselves and arbitrary 
other objects. The pay for this simphfication is that DS evolution as determined by Eqs. 16 
and 18 becomes non- unitary, because of r]{t) 7^ C(^) ■ It is useful to introduce new noises, 

^ = ,x + iy/2 , r] = x- iy/2 , x = {7] + 0/2 , y = i{r] - ^ 

In their terms Eqs. 17 and 18 read 

/exp fj2[9jit)xjit) + Mt)yjit)]dt) = (19) 



9 At) 



+ ifj{t) 



Bj{t)dt Sep / ^ 



9 jit) 



-if jit) 



Bjit)dt } pb 



with gj{t) and fj{t) being test functions, and Eq.l8 takes the form 
d 



dt 



p{t)^v{t)Jop + L{t)p{t), L{t)p^J2yjit)Djop + i[p,H{x,t)] , 



(20) 



where H{x,t) is defined again by replacing B{t) in (11), H{x,t) = Heit)+ J2j Xj{t)Dj . 
According to (19), in bath phase space x{t) is associated with time-ordered symmetric mul- 
tiplying, B{t) o ... , while y{t) with time-ordered commutation, i[B{t), ...] . 

At v{t) — , Eq.20 (or equivalently Eq.l8) turns into what can be named generahzed 
Langevin stochastic equation for DS interacting with a thermostat. It differs from frequently 
used quantum Langevin equations by the terms containing yjit) . Clearly, just these terms 
destroy unitarity of marginal DS evolution. Like the testing they represent non-Hamiltonian 
perturbation of DS. Their physical meaning is the phase volume exchange between DS and 
bath. Of course, the unitarity effectively restores after the complete averaging is performed. 

Notice that evidently dissipative terms appear in Eq.20 if divide products Xj{t)p into the 
averages {xj{t)p)^ and normalized forms : Xj{t)p : . From the other hand, comparison of 
Eqs. 10 and 20 shows that Ti gp{t) may be interpreted as mutual characteristic functional 
of observables J{t) and Dj{t) as if governed by Hamiltonian H{x,t) and tested by v{t) and 
yj{t) , respectively. 

It should be emphasized that Eq.20 also directly extends to classical mechanics if replace 
commutators by classical Poisson brackets. From the point of view of classical intrinsic bath 
dynamics, multiplying by y{t) looks as chronologically ordered action of the brackets. 



yit) ^ E 



dB{t,q,p) d 
dp dq 



E 



dB{t,q,p) d 
dq dp ' 



where q,p are bath canonical coordinates and momenta, B{t, q,p) are B{t) as expressed via 
initial values of q,p , and differentiations d/dq , d/dp are oriented in the direction of initial 
bath distribution pb — Pbiq,p) ■ 
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G. Effective bath noise statistics 



Thus yj{t) substitute superoperators even in classical case. Nevertheless, even in quantum 
case we may treat them as well as x{t) like c- number variables if only remember the definition 
(19) and carefully keep in mind their time assigning. The evident profit of such a symbolic 
calculus is that now a "ready-made" statistics of bath noises x{t) and y{t) may serve as a 
model fit of the interaction. But underlying operator nature of y{t) results in their unusual 
partly "ghost" effective statistical properties. For instance, at a{t) = Eq.l9 turns into 



ewljY^bj{t)yj{t)dt\\ =1, 



that is all statistical moments of yj(t) are zero. In spite of this, yj{t) are non-zero, since they 
have non-zero cross correlations with Xj{t) . 

Of course, x{t) and y{t) are closely related being two manifestations of the same intrin- 
sic bath life. General connections between them follow from the properties of the "trace" 
procedure. In the stationary case, when [Hi,, ph] = , for two-point correlators (which are 
translationally invariant in this case as well as higher correlators) we have 

Kj^ir) = IV, S,(t)S^(0) p, = K*^j{-r) (21) 

and correspondingly from Eqs.l7 and 19 obtain 



(xx) 



T = [Xi 



T)xmm = ReK^rnir) ; i^j^ (r) = (a;,(r)2/„(0)) 



2 Im Kjm{r) , r > 
, T < 



(22) 



Here asterisk denotes complex conjugate, and subscript "6" is omitted. 
Under the canonical bath distribution, 

Pb oc exp(-if(,/9) , 

where © is bath temperature, additional relation takes place, 

Kj^{T-i/e) = Kmj{-r) (23) 

The latter if combined with (21) means that 

-^im V) = cos{uJT)ajmi(jj)duJ , i^j^^V) = 2^?('^) sm{ujT) tanh (^^^ aj^{(jj)duj 

with some real non-negative spectral measure ajm{u!) and ■i?(r) being step function. 

In the important specific case, when B(t) represent a free bosonic field and therefore 
self-commutator of B{t) is c-number. 



7 



i[Bj{t), Bm{t')] = 2Im Kmj{t' - t) , 
for any functional ^{x,y) the Eq.l9 implies 

{yjim^, y)) = ( / E[2 Im K^j{t' - t)] 

\Jt'>t ^ 

In this sense, in case of free bosonic thermostat 



(24) 



H. Ghost fields and dissipation 

According to relations (22) and (24), all ghost noises y{t) are correlated with more late 
x{t) only. The meaning of this property becomes trivial if rewrite Eq.l9 in the form 



lY^[9j{t)xj{t) + fj{t)yMdt) = 



where Bj{t, f) correspond to non-free bath evolution perturbed by classical forces f{t) , 



Hb{t')dt' 



Bj^xp 



-I I Hb(t')dt' 





, H,{t)^H,+j2fjmj 



Differentiation of this identity by g{t) and f{t) shows that cross correlators are nothing but 
bath response functions to the perturbation: 



n 
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\ j m I 

Consequently, for all higher-order correlators 

(n^fe) n^(4) ) = if max{4} > max{tj} , 



(25) 



/=o 



(26) 



in accordance with causality principle. This property just guarantees eventual safety of the 
unitarity and of summary phase volume of two subsystems. Due to (26), as one can see from 
Eq.20, at v{t) = always Tre {p{t))^ = 1 . 

Because of obvious dual symmetry of the subsystems, with reference to bath f{t) and g{t) 
play in Eq.25 the same roles of exciting forces and test functions, respectively, as x{t) and 
y{t) in Eq.20 with reference to DS. It is not hard to guess that the intention of variables y{t) 
in the symbolic bath noise calculus is involving a back action of DS onto bath, i.e. reverse 
dissipative energy flow which must ensure conservation of total energy. If y{t) were absent 
then only non-reciprocal pumping of DS by forces x{t) would be present. 
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I. Fermionic dynamical subsystem 



Let DS is purely fermionic representing free conducting electrons (plus some electron 
depository if necessary). Our main ansatz will be that (vi) all the operators related to this 
subsystem have the quadratic form 

A = Y.A.pCiC(, (27) 

where (C) are fermionic creation (destruction) operators, indices «, (3 enumerate some 
complete set of one-fermion states, and A stands for i/g , Dj , Q and consequently J . 
Thus the Hamiltonian (2) does not contain direct electron-electron interactions like the 
Coulombian one. But this simplification implies no principal losses, since it is known that 
such interactions can be introduced by means of auxiliary actual or ghost bosonic degrees of 
freedom included to both Hb and the set Bj . 

Indeed, let Hg is quadratic and we add to it " Coulombian" or any other pair interaction 
described by the Hamiltonian 

ad 

(of course, here a, j3 may relate to different set of electron states than in (27)). Then Eq.20 
turns into 

j^p{t) = v{t)J o p{t) + L{t)p{t) + i\p{t), H,] 
It is easy verifiable that its solution can be represented in the form 

pit) = m) , = v{t)J o p{t) + L{t)p{t) + L\t)p{t) , (28) 



L'{t)p = Y.y'ama o-p-riY^x'^m, , 

a a 

where x'{t) and y'{t) are Gaussian white noises with correlators 

{x\t)x\to)) = {y'{t)y'{to)) = , {x'^{t)y'^{to)) = 2VapS{t - h) , (29) 

and angle brackets in (28) notate averaging with respect to these white noises if mentioned 
in the sense of Stratonovich (not of Ito). 

This is merely the variant of so called Stratonovich transformation. Now all fermionic 
operators in the evolution equation (28) become quadratic. Hence, after including x'{t) and 
y'{t) to the set of bath noises (with operators Na in place of Dj ) electron-electron interaction 
reduces to the above general stochastic scheme. Because of {x'x') — , one may say that 
x'{t) and y'{t) represent interaction with zero-energy bath. 
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J. Fermions to matrices conversion 



To complete our simplifications, in addition to (14) we assume (vii) an initial distribution 
of the electron DS be canonically produced by quadratic-type operators. For instance, 

^ exp(-go/T) exp(/x Eg N^/T) 

IVe{exp(-//o/T)exp(/.EaA^a/T)} ' ^'"'> 

where Hq is some quadratic operator, and sum X^^, is the operator of total number of 
electrons (it commutes with any other operators). In particular case, Hq may coinside with 
Hf, . If wanting to fix a number of electrons, Nf.i , instead of chemical potential, ji , we may 
consider the composition of grand canonical distributions, 

Tr,{<5(iVe,-E«iV„)exp(-i7o/T)} ' '^""^ 1/ 27r ' ^"'^ 

with variety of imaginary chemical ponentials, itpT . In this approach, thermodynamical 
effects of Coulombian interaction, if included into dynamics as prescripted by Eqs.28 and 29 
(but with the averaging put off to a later stage), will come into play after a time of evolution. 

Then, clearly, in Eq.l6 we can confine ourselves by calculation of only specific fcrmionic 
spurs Trg {Hn^n} I where each multiplier is an exponent of one of three sorts, exp(yl) , 
^xp{/ A{t)dt} , or exj3{/ A{t)dt} , with A being quadratic fermionic operator in the sense of 
(27). This calculation is easy performable if take into account three remarkable facts, namely, 

(1) any exponent of the second or third sort can be written as an exponent of the first sort, 

(2) product of two such "quadratic exponents" can be represented as a single quadratic 
exponent and (3) there is one-to-one correspondence between such operator expressions and 
analogous matrix expressions. In summary, 

exp(^X,^C+C^)exp(^r,^C+C^)=exp(^Z,^C+C^) ^ e^e^ = (32) 



6xp 



j Y.X^p{t)C+C0dt^ = exp(^Z„^C+C^) ^ gxp [1 X{t)dt 



(33) 



Here Xaii = (^)a/3 , and so on, and the right-hand sides are meant in the sense of matrix 
functions. These relations can be proved if use the commutator identity 

Hence, any product of quadratic fermionic exponents can be transformed into a single 
such exponent, for which it is known that 

TVe exp(E Z^pC+Cfs) = det (1 + e^) = exp Sp ln(l + e^) (34) 

with det (...) being matrix determinant, and Sp M = J2a ^aa for any matrix M . After that, 

according to (32) and (33), the exponent may be back unfolded as the exact matrix copy 
of initial operator product. Therefore, if denote quadratic operators and their exponents by 
the same letters as their matrix copies then 

Tre {H Mn} = det {1 + n (35) 

n n 

Although x{t) and especially y(t) are to some extent ghost variables, we can pass them 
through all these transformations like scalar c-number variables, since they by their definition 
commute with arbitrary objects and are uniquely marked by their time arguments. 
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K. Secondary randomization instead of secondary quantization 



Let po be matrix copy of quadratic fermionic exponent ex]){—Ho/T) , with Hq now defined 
in the space of one-electron states. Then, as it follows from previous section, in case of grand 
canonical ensemble (30) the characteristic function (16) can be represented by 



(det[l + zpjt)]), _ (cxp Sp ln[l + zp{t)])^ 
det(l + zpo) exp Sp ln(l + zpo) 



z = e^l^ 



where time-dependent one-particle density matrix 



pit) = gxp 



-v{t')J-tH{^,t') 



dt'^ Po ex^ i^j^ 



-v{t')J + iH{r^,t') 



dt' 



(36) 



(37) 



is matrix copy of the above considered p{t) . It obeys exact matrix analogies of Langevin 
equations (18) and (20), but now under initial condition p(0) = po . The case of canonical 
ensemble (31) can be reduced to calculating the same quantities as in (36): 



E{v) 




(exp Sp ln{l + zp{t)}) I 



J z=0 L 




exp Sp ln(l -|- zpo) 



(38) 



2=0 



Now, when Fermi statistics is explicitly took into account, we may pay our attention to 
one-electron version of Langevin equations (18) and (20) and analyse statistical properties 
of one-electron density matrix (37) (or corresponding wave functions) as evolving under 
influence by random processes x{t) and y{t) . If we first performed bath averaging we would 
deal with evolution of many secondary quantized fermions. Instead we will deal with twice 
random (quantumly and excitedly) motion of single particle. 

The above consideration clearly demonstrated that we have all rights to deal with x{t) and 
y{t) as real- valued variables and correspondingly with 77 (t) and ^{t) as complex conjugated 
variables, ri{t) = (C(^))* • Hence, if suppose test function v{t) be also real- valued then we 
can formally treat p{t) as Hermitian matrix and write 



p{t) ^ S{t,to)p{to)S+{t,to) , S{t,to) = ^(^j' ^v{t')J -iH{C,t' 



dt' 



(39) 



In general, S{t, to) is not unitary transformation, SS~^ ^ 1 . 



L. Grand canonical ensemble 

Consider general relations what follow from one-particle matrix versions of Langevin 
equations (18) and (20). It is easy to see from Eq.20 that 

for any function F{p) — F{p{t)) . Therefore, if putting on F{p) — ln(l -|- zp) , we may 
rewrite Eq.36 as 
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E{v) = (exp 



I v{t) Sp J fit) dt + Y.j vAt) ^vDjfit) dt 



, fit) 



1 + zp{t) 



(40) 



Clearly, matrix f{t) reflects occupancy of one-electron states in factual many-electron system. 
According to Eq.20, fit) is solution to the non-linear stochastic equation 

f = ^[/^(l -/) + (!- f)Jf] + E ^-^[fm -/) + (!- mf] + t)] (41) 

To find statistical moments of the current we have to differentiate '^iv) by test function 
v(t) and then turn it into zero. But the terms in (40) containing fjjit) remain and work even 
at u(t) = . Hence, it is comfortable to introduce special notation for bath averaging under 
vit) = : 



(*(x,y))o = 



^ix,y) exp ^ / yjit)SpDjfit)dt 



v=0' b 



with $(a;, y) being arbitrary functional of bath influence. Due to the property (26), (1)q = 1 . 
Because of (26), here as well as in (40) in the exponent we may omit limits of the integration. 
Then the mean current is 



{Jit))^mv)/sv{t)u,^{spjfit)), 



(42) 



The current correlation function produced by second-order derivative of S(v) generally looks 
more complicated: 



(J(t)J(to)) = S^Eiv)/dvit)5vito] 



v=0 



(43) 



Sp {fit) J} Sp Sit, to) [J o zpito)]S+it, to){l + zpit)}-')^ + 



-1 



[SpJ{l + zpit)}-^Sit,to)[Jozpito)]S+it,to){l + zpit)y , 
where all matrices correspond to vit) — , and it is assumed that t > to . 



M. Non-degenerated electrons 

Let number of electrons is small as compared with total number of one-electron states, 
Ngi << N , N = Sp 1 , to the extent that all mean occupancies are significantly smaller 
than unit. Then, in the frame of canonical ensemble, Eq.38 approximately transforms into 

Eiv) ^ ([Sp pit)r-)^ , Pit) ^ (44) 

Here the normalized density matrix pit) , of course, also obeys Eqs.18 and 20. In its tems. 
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(Jit)) = Nei (Sp Jp{t))^^^_, (45) 
{Jit)J{to)) = Nei{N,i-l)(Sp{Jp{t)}Sp{S{t,to)[Jop{to)]S^it,to)})^^ (46) 
+Nei (Sp JS{t, to)[J o p(io)]^+(i, to)] 



Nei-2 



where t > to , and now the bath averaging at v{t) =0 is defined by 

mx,y))^ = ($(x,|/) [Sp p{t)r), ^ U{x,y) expnJ2j vAt) SpD,p{t) dt\ 

\ J ' b 

(the latter expression relates to large N and n under small ratio n/N « 1 ). 

Since generally Sp p{t) 7^ 1 , we see that even in dilute electron gas contribution of 
any separate electron to total current and its fluctuations statistically depends on all other 
electrons. The matter is that all they feel the same thermostat and thus the latter becomes 
mediator of electron-electron correlations if not interactions. Since non-degenerated electrons 
obey Boltzmannian statistics, formulas (44)-(46) relate also to classical hmit. 



N. Diagonalization of stochastic density matrix 

According to Eqs.36, 38 and 44, in principle, to find the characteristic functional we need 
in eigenvalues, {k = 1..N), of the one-electron density matrix only. In view of Eq.39, we 
can treat p(t) as Hcrmitian operator and write p(t) = \E'(t)A(t)\E'"'"(t) , where A = diag{Afc} 
with real elements and "^{t) is unitary transformation. Then in grand canonical ensemble 

llfe=i[l + ^AofeJ 

where Wx{t, A) = Wx{t, Ai, .., Ajv) is probability density distribution of the eigenvalues, dX — 
rijfcLi dXk , and Xok are their initial values. In ensemble with a fixed number of particles 

"(^) = / J2 Afe,Ajk2..Afc^^ 1 iyA(^,A)dA AofeiAofe2-Aofe;v,J ~ (48) 

\ki<k2<- ^ I Vfei<fe2<... " I 



~ / E-^M ^A(t,A)rfA^Aofc 
^ \k=i / \fe=i 

where the latter expression relates to dilute gas. 

To realize this most adequate approach to the problem, we should combine evolution 
equation (20) and the first-order perturbation theory to describe infinitesimal time steps of 
the diagonalization. After standard manipulations we obtain coupled evolution equations 
for eigenvalues and eigenfunctions of density matrix: 
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dt 



it)i^l,j^k) + Eyjii)in,D,^k) 



(49) 



d-^k 
dt 



it)i^*^,J^,) + Y.y,it)i^*^,D,^,) 



-iH{x,t)^k (50) 



Here (a, 6) = Z^Oa^'a means pseudo-scalar product of two complex vectors (thus (a*, 6) is 
usual scalar product) , and we suppose that eigenvalues are non-degenerated. The latter 
assumption is not too restrictive since will acquire tendency to mutual repulsion. 



O. Weak noise and Fokker-Planck equation 

Let us emphasize that in fact the Langevin stochastic equations (49) and (50), as well 
as Eqs.18 and 20, are still exact equations, in the sense that they represent explicit mould 
of underlying detail Hamiltonian dynamics of the electron-bath interaction. The natural 
next stage of the investigation is derivation of corresponding kinetic equation for probability 
distribution of A and \I' (and thus for probability functional of p{t) ). 

If x{t) and y{t) were white noises the result would be time-local Fokker-Planck (or Kol- 
mogorov) equation again exactly reflecting the Hamiltonian dynamics. As Eqs.28 and 29 
demonstrate, this is just the case when the noises represent time-local electron-electron inter- 
action only. We put off this interesting example for the future. Another such case is y(t) = 
while x{t) being white noise which may be interpreted as the infinite bath temperature limit. 
We will touch it below. 

In general, exact kinetic equation is inevitably non-local since at any finite bath temper- 
ature the relation (23) forbids to at once turn all the correlators (22) into delta-functions. 
But we believe that a time-local kinetic equation will be as usually good approximation if 
the bath noise is sufficiently weak and short correlated. Moreover, under such assumptions 
we can use simplest one-loop approximation resulting in a Fokker-Planck equation. For any 
set of stochastic equations 

dZ 

^ Ak{t,Z) + Y,Xa{t)Bak{Z) + ^ya{t)Cak{Z) , 

with the same structure as Eqs.49 and 50, this approximation yields 

dW d d 

= Ait)W = Afit)W - g^F,{t) W , Af{t) ^-Y. g^A,{t, Z) , 

where W = W{t,Z) is probability density function, operator Af{t) corresponds to free 
evolution and operators F^it) describe the bath noise contribution to the probability flow, 

poo /f 

F,{t)^-YBakiZ) / (x.(r)x^(0))e-^/W-^S^„(Z)e--^^Wdr- 

Jo OZ/m. 



- Y B^,{Z) r (x„(r)y^(0)) e-^/W-|-C^„(Z)e-^^/(*)dT 
(supposing noises are stationary and have zero mean values). 
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p. Kinetics of probability measure of density matrix 



When applying this scheme it to our Eqs.49 and 50 we may neglect a contribution by v{t) 
to the exponents like exp[±r A/(t)] in the above integrals, since under weak and short corre- 
lated bath noise this contribution might result in relatively small corrections only (physically, 
this means that we neglect displacements of an electron during its short collisions with bath 

excitations). 

Then Fokker-Planck equation for the mutual probability density distribution of the eigen- 
values and eigenvectors, W{t, A, ^) , looks as 

dW 

= K{t)W = [v{t)Kt,st + A/ree(^) + A..(^) + Ky{t)]W , (51) 

where kinetic operators latest , ^free , ^xx and A^^y describe testing of electron DS, its free 
evolution as if in absense of noise, its stochastic pumping by x{t) which itself influences 
only, and common action of y{t) and x(t) which leads to dissipation and besides to 
fluctuations of , respectively. These four parts of A{t) directly correspond to four different 
type terms in stochastic equations (49)- (50). The first is 

' ' d 



(52) 



with d/d'^k being gradients and {d/d'^k)* = d/d^l . 

To express other parts of the kinetic equation (51), for any Hermitian matrix M let us 
define first-order differential operator 



n{M}^zE 



d 



d 



Obviously, Afree{t) is such operator, Afree{t) = ^{He{t)} , while effect of noise x{t) is 
described by the second-order operator 



Kx{t) = E^{^a}^{ / K^:;\r)exp[-iH,{t)T]D^eM'^H,{t)T]dT 

The second-order operator Axy{t) contains x-y-correlators instead of x-x ones, 

d 



(53) 



A^y{t) = -Q{L»J^A, (vl>*, J^^ (r)e-^«W^i?^e^^«W^dr*,) - (54) 



d 



with summation over all repeated indices and k ^ m . Its structure resembles that of Atest 
and it also mixes kinetics of ^ and A . 

What is for the case of sufficiently strong DS-bath coupling (intensive bath noise), one 
can make the substitution 

ft 



H{t',x)dt' 



in Eqs. 49-50, which allows to treat (^J^, S^{t)DjS3;{t)^k) and similar quantities as fast noises 
and again apply Fokker- Plank ideology. 
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Q. Informational vacuum 



Let Wunifi"^) be uniform distribution of ^ which is invariant with respect to arbitrary 
rotations and thus independent on ^ at alL Of course, this uniform measure is concentrated 
on the unit sphere, {^^, ^k) — ^mk ■, as well as W{t, A, ^) in general, as Eq.50 does imply. 

It is not hard to verify two facts. First, multiplication by Wunif commutes with complete 
kinetic operator A(i) and with any of its four parts separately, 

KWunifi^) - Wunif{^)K = , s e {test] free] xx] xy} (55) 
Second, kinetic equation (51) always has formal stationary solution 

Wstat{\ *) = W^^if{-^)Y{\) , IV{t)Wstat{\ *) = , (56) 

where function Y{\) is common solution to the equations 

+ E yi>^) = , = (57) 

The latter reads 

(Afc — A„) 

llj=l 

Moreover, each of the parts of A(t) separately turns Wstat{\ ^) into zero, 

KWstat{\ *) = , s e {test] free] xx] xy} 

The immovable point Wstat{\ ^) can be named informational vacuum since it says nothing 
about the system. 



R. Much more stationary states 

However, this stationary solution can be attracting point for solution of Eq.51 at non-zero 
test parameter, v{t) ^ Q , i.e. in absence of testing only. To see this let us integrate both 
sides of Eq.51 over ^ . Since the only part of A(t) what has terms free of gradients d/d'^k 
and d/d^% is A^st , only its contribution survives after this integration: 

-Wx{t, A) = -v{t) Y: ^Afe I {%, J^k)W{t, A, *)d* , (59) 



Wx{t,X) = j W{t,X,^!)d^! 

Here Wxit^X) is marginal probability distribution of eigenvalues. At v{t) = it does not 
change at all and remains atomic, W\{t,\)\v=Q = S{X — Aq) , with Aq being initial values. 
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Consequently, under a constant perturbation, u{t) — const , W{t, A, tends to a stationary 
distribution, 

iyoo(A,*|Ao) = lim W{t,X,^)l=o , 

t—*oo 

which differs from (56) and remembers exact start eigenvalues. Hence, at v{t) — there is 
continuous infinite set of stationary solutions to Eq.51. 

Interestingly, in general these degenerated solutions are not atomic with respect to A , 
i.e. H^oo(A, ^'|Ao) 7^ S{X — Aq) • J Woo{X,'^\\o)d'^ . This unusual property as well as the 
degeneracy itself is reflection of ghost properties of y{t) . The latters manifest themselves 
in absence of second derivatives d'^/dX'^ in Eq.51, in spite of presence of joint derivatives 
d"^/ dXd^ , 9^/ dXd^* . Therefore, strictly speaking, A{t) is not purely elliptic operator. 



S. Exclusion of eigenvectors 

In order to analyse the characteristic functional as presented by Eqs.47 or 48, we should 
consider marginal distribution Wx{t,X) just ai v{t) ^ , when it can change to arbitrary 
extent. Since then the vacuum solution (56) is attractive, it is natural to use it as ground state 
in the projection method. All the more, bath noise also pushes W {t, A, ^) towards Wstat{X, 
making eigenvectors as much random and statistically independent on eigenvalues as possible. 
At the same time, since the factor Y{X) is non-normalizable, this factored solution never can 
be achieved explicitly. Always more or less correlations between A and do exist, and we 
must write 

Wit, A, ^) = W^it, X)Wunifi^) + Wit, A, *) , / Wit, X, ^)d^ = , 

where Wit, X, "if) does not contribute to Wxit, A) . 
Thus we merely introduce the projection operator: 

UWiX, *) = W^r^^fi^) J W iX, ^)d^ , Wit, X, ^) = (1 - H) W(t, A, *) 

With using Eqs. 55-56, Eq.51 yields 
d 



mWx- I WunifMt)Wxd^ 

Here factor /(t) in the source term is easily calculable: 



gW^il- U)Ait)W + Wuniflit) , lit) = 



(60) 



d ^ . sr^ Xm + Xk 



m = - E + E Wxit, X)[vit)i%, J*,) + i%, F(t)*,)] , (61) 

with the traceless matrix Vit) defined by the sum of commutators 

Vit) ^iY.\r K^:^\T)eM-iHeit)T}DpeMiH,it)T}dT ,D^ (62) 

Only Ktest and K^y contribute to the source giving two similar terms in (61). Wee also took 
into account that 



j (**,M*fe)W^„,,^(*)ci* = (SpM)/Ar 



for any matrix M . Because of this equality and of that the current matrix J by its definition 
is traceless, only W contributes to right-hand side in Eq.59. 
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T. Kinetic equation for eigenvectors 



As usually in projection techniques, Eqs. 59-62 ensure exclusion of a part of variables, 
producing formally closed although non-local separate evolution equation for eigenvalues. If 
introduce notations 

A,, ^ {n,A^,) , ((...))_^ ^ / {..^w^mmd"^ , 

Uvit.to) = Sep [V{t'){l - U)Atest + Afree{t') + Axx(t') + A^y{t')] dt'j , (63) 

this equation reads 

^W,{t, A) = v{t) Y: ^Afe 1^ {hkfJvit, t') [v{t')Jqq + yM)^^,f X (64) 

X (^^«+ E I W^{t',\)dt'-v{t)Y.^\, I J,fcC/,(i,0)ty(0,A,*)d* 

We take in mind Eqs. 55-56 and the identities (1 — n)As = A^ for s e { free; xx; xy} . 
The important point is that last term which represents a legacy of initial state, 

W{0, A, ^) = 6{X - Xo)m - ^o) - Wum/m] , 

may be insignificant even in spite of the degeneracy. The matter is that all the appointed 
stationary distributions do not depend on initial ^ distribution, therefore Uy(t. 0) W (0, A, ^') 
must decay after a fixed finite time determined mainly by the diffusion operator Aj.x , i-e. 
characteristic correlation time of eigenvectors fluctuations. At the same time, as already 
Eq.59 shows, the evolution rate of Wx{t, A) is dictated by magnitude of test function v{t) . 
Since the latter serves eventually for differentiation by it, we may treat it as infinitesimally 
small parameter and hence Wx{t, A) as slow varying, thus having all rights to neglect the 
last term. 



U. Kinetics of eigenvalues in hot bath 

Let the perturbation u{t) which namely determines time dependence of operators A free 
, A^x ^-nd Axy is not too strong and not too fast varying (for instance, u{t) is constant). 
Then we can neglect time variations of u{t) in the integrand kernel in Eq.64 while in (62) 
neglect u{t) at all. Besides, let us neglect also the kernel dependencies on both v{t) and A , 
i.e throw away Atest and A^y from the propagator (63). Formally, it is possible because the 
diffusion operator A^x itself quite ensures decay of the kernel, while effects of operator A^y 
is partially took into account by the term containing (^^,y(i)^fe) in the source (61). This 
term involves dissipation and represents the storage for non-zero dissipative mean current 
response. 
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Clearly, in such an approximation the current is modelled as Gaussian random process, 
and all effects of the degeneracy are asided. Neglecting K^y as compared with A^^ corresponds 
to what can be named "infinitely hot bath limit", since according to Eq.23 action of K^x 
much overpowers action of h.xy when bath temperature © grows to infinity. Thus we come 
to maximally simplified kinetic equation for eigenvalues distribution, in the form 

^^^?^-(^-^?.^)- 

Here the evolution rate V{t) is determined by 

m 



N 



poo 

,it) I {i%,J^o)eM[^freeit)+AMi%,V^o))^^,fdT+ (66) 



roo 

+V{t) J^o) exp{[Ay,eeW + Axx]r}{%, J^o))unif<^ ' 0^?^ ' 

with Kfree{t) = ^{Hg} — u{t)Q{Q} , vector is any of columns of , and matrix V is 
given by (62) but under He{t) replaced by He . We assumed also that total number of 
one-electron states (dimensionality of matrices) is large, N » 1 , which allows for many 
formal simplifications. 

The first integral on right hand of Eq.66 describes mean current. Wc may suppposc 
that this integral turns into zero at M(t) = , because of a proper symmetry of the matrix 
(62), i.e. that there is no spontaneous current breaking time symmetry. We expect also at 
non-zero low-field conductivity and sufficiently small and slow varying u{t) (and besides slow 
v{t) ) expression (66) reduces to 

T{t) a; v{t)^u{t) + Av'^{t) = v{t)w{t) + Av'^{t) , (67) 

where A and w{t) represent diffusivity and drift velocity of a separate electron. 

Notice that Eq.65 closely resembles kinetic equations for transmission matrix eigenvalues 
in the theory of statically disordered conductors [16]. 



V. Current noise in Gaussian approximation 

Rather trivial model (65)- (67) is nevertheless useful probe to compare roles of diffusive 
component of Wx{t,X) evolution (described by second-order derivatives d'^/dX^ in Eq.65), 
from one hand, and of drift caused by mutual repulsion of eigenvalues (described by ratios 
i^m + ^k) /i^rn — ^k) )• from another hand. Note that at possibly negative r(t) values Eq.65 
corresponds to mathematically incorrect time-reversed diffusion, but in fact this implies 
nothing incorrect in final results of calculations. 

Simple analysis shows that, at >> 1 , from the point of view of statistical moments of 
A , all diffusive terms in Eq.65 occur relatively insignificant. Therefore, in the asymptotical 
sense, for any function F{X) we can write 
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It 



W)) = r(t)(^E|E 



Afc + A, 

Afc — At 



;) A. A^(A)), ((...)> ^/(...W 



(68) 



This circumstance is due to that many of ratios (A^ + A^) / (A^ — A^) are large numbers 
~ N , hence the drift overpowers diffusion approximately N times. 

Particularly, in case (48), i.e. at fixed N^i (may be comparable with N ), 



EE 



Afc + A^ 

Afc — A,; 



d ' ' 
'^'^^^ E '^'=i'^'=2--'^ifcjv^, — ^ei{N — Ngi) ^ AjfcjAjk2--Afcj^^^ (69) 



(primed sum means that all indices are different one from another), and Eq.68 exactly yields 



E{v)^ew[eNei{N-N,i)/N] , 6= f V{t')dt' 

Jo 



(70) 



This result corresponds to Gaussian current noise, with only those electron-electron statis- 
tical correlations what come from Pauli principle and Fermi statistics of electrons (or holes 

In case of grand canonical ensemble (47), we obtain 



exp 



N^N-Nei) , 9 



N 



+ 



1 - 



N 



EA(i-^) + 

k 



(71) 



fk = zXk/{l + zXk) , Nel = Y,fk, 

k 

where A^ should be replaced by their initial values (i.e. represent diagonal elements of 
p{t = 0) ). The points mean higher degrees of 9 . The second term in the exponent 
has no analogy in (70) and formally corresponds to non-decaying flicker type current and 
conductance fluctuations. But in fact this is nothing but artifact of grand canonical ensemble. 



W. Degeneracy and low-frequency noise 

Now let us be so courageous and discuss what will occur beyond the hot bath Gaussian 
approximation, i.e. if keep v , A and d/d\ dependencies of the kernel in Eq.64. Generally 
speaking, we may expect that above mentioned strong degeneracy of stationary state results 
in long-lasting statistical correlations when observing DS, i.e. in low-frequency excess fluc- 
tuations of the current. Indeed, as Eq.59 shows, measuring of the current transforms initial 
values Ao to one or another different values A' which then remain for arbitrary long time 
being unchangably encoded in iyoo(A,^|A') (or in similar time varying distribution if u{t) 
is not constant). Therefore, next measurement will be inevitably correlated with previous 
one regardless of their time separation (this resembles low-frequency Goldstonian excitations 
under breaking a continuous degeneracy). 

If this reasonings are true, then the differential operator kernel in Eq.64, 

Rkg{t, to, V) = V{t) {jkkU,{t, to) [V{t0)jgg + Vgg{to)]) , (72) 
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has slow decaying non-exponential tail closely related to continuous degeneracy of ground 
state of the operator (1 — n)A(t) which generates semi-group U^,{t, to) . Clearly, such a degen- 
eracy takes place (and besides, in opposite to A{t) , even at arbitrary v{t) ), hence, Uy{t,to) 
is not strictly contracting semigroup. Directly, the degeneracy concerns A dependence of 
joint probability distribution, but indirectly it expands also to its ^ dependence, because 
of the \ - coupling produced by operator A^yit) . The hot bath Gaussian approximation 
considered neglects this coupling and therefore retains only rapidly decaying part (66) of 
the kernel (72) responsible for mean diffusivity and mean drift. The rest of the kernel must 
contain information on excess contribution to the current correlation function, 

K{t,to) = {J{t)J{to))-{J{t)){J{to)) , 

including conductance fluctuations, and besides on higher-order cumulants of current. The 
former, in its turn, consists of two parts whose physical sense is as follows. 

When the bath noise x{t) transforms coherent electron motion into diffusion (plus drift as 
characteristics of its asymmetry), diffusivity and mobility of any electron depend on random 
time-varying non-uniformities of density of other electrons, i.e. on random occupancies of 
one-electron states. As it will be demonstrated elsewhere, this fact results in flicker type 
conductance fluctuations with spectrum oc f~^^^ , at frequencies higher then inverse time 
of electron diffusion through overall system. This excess noise arises merely from Pauli 
principle and has no relation to A - ^ coupling. The only role of Axy{t) in this mechanism 
is establishing more or less statistical hierarchy (non-equiparticity) of electron energies, in 
contrary to x{t) pumping, and thus translate diffusivity fluctuations into drift ones. 

Other mechanism of excess low-frequency noise is due, in principle, just to the A - ^ 
coupling. As we already pointed out, the reverse DS action onto the bath, representing 
dissipation, makes the latter mediator of indirect electron-electron interaction through bath. 
This is in rich analogy with action of the x'{t) - y'{t) noise in Eqs. 28-29 which is equivalent 
to direct electron-electron interaction. Hence, to keep operator Aa,y(t) in (72) means to take 
into account such an interaction. But our experience gave evidents [5,8] that particle-particle 
interaction itself is sufficient natural cause of a true flicker noise, i.e. noise without low- 
frequency spectrum saturation (1/f fluctuations of self- diffusivity and mobility of particles). 
In [1,2,4,5] it was explained why this noise can be insensible to spatial and temporal scales 
of a system. 



X. Current correlation function 

To connect K(t,to) and the kernel (72), let us suppose, for simplicity and visuality, that 
external perturbation u{t) is time-independent and that after not a long time a mean value 
of the current, J = {J{t)) , tends to a constant, while its correlation function becomes 
depending on time difference only, K{t — to) . Then we can put on the test function v{t) 
also constant, v — const , thus testing integrated current, i.e. transported charge, /q J{t')dt' 
. This simpliflcations bring comforts of Laplas transformation. 

Applying the latter to the characteristic functional, at sufficiently small v and p we have 

E(v) = / E(v)e-P'dt = \p-vJ - v^Kip) - ...]-^ , Kip) = / KMe-^^dr , (73) 
Jo Jo 
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where the dots replace higher-order terms of v expansion. From the other hand, the same 
quantity can be expressed from Eqs.47-48 and evolution equation (64) (throwing out last 
term of (64)), 

^v) = F-\X) [p - vGiip) - v'^G^ip) - ...] ' F{\) (74) 

Here differential operators Gnip) acting on A are defined by v expansion of kernel (72), and 
A stand for the initial eigenvalues. Direct formal expansion gives 

Gi{p) = YXd, + '/'^(A)] {Jkk\p - Ao , 

kq 

G2{p) = Y^idq + 0,(A)] Uk\p - Ao ]~'^) dk+ (75) 

kq 

+ T.[9, + M^)] {Jkklp - A^]-^(i - n)ALb - A^]- V,,) d, , 

kq 

where the shortened notations are used 

dk = ''^k-^T- , 0fc(A) = ^ ' Ao = A/ree + ^xx + ^xy 

The superscript "T" is symbol of conjugation with respect to A in the Sturm-Liouville sense, 
that is reversing of order of all differentiations d/dX and multiplications by any function of 
A and changing d/dX by —d/dX . 

Comparing linear and quadratic terms in formulas (73) and (74) we see that 

J ^ F-\X)G,{p)F{X) , (76) 

K{p) = F-\X)G,ip)FiX) +p-'{F-\X)Gi{p)F{X) - [F-^(A)Gi(p)F(A)]2} 

Thus, mean current really tends to a constant if Gi{p 0) has finite limit. Besides, the 
natural possibility is if in this long-time limit F{X) is eigenfunction of operator Gi{p 0) . 
Then Laplas transform of the current correlator reduces to 

K(p) ^ F-\X)G,(p)F(X) (77) 

Hence, F{X) may be also eigenfunction of G2{p 0) . Notice that operators (75) always 
can have polynomial eigenfunctions like presented by (48) since these operators (as well 
as complete operators Aq , A^gst ) are insensible to arbitrary scaling of A . Of course, 
both J and excess part of K{p) are results of external perturbation, u , hidden in A free — 
n{He} - un{Q} . 
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Y. Resume 



For the meanwhile, our main results are, first, reduction of many-electron problems to 
investigation of stochastic one-electron density matrix, p{t) , which obeys exact Langevin 
equation (20), and, second, formulation of approximate Fokker-Plank kinetics for eigenvalues 
and eigenvectors of this matrix. 

Equivalently, one could consider statistical moments {p(t) ® ... ® p(i))f, , i-e. mean strict 
products of many copies of p(t) , and derive from (20) evolution equations for them or, 
instead, express solution to (20) and then the moments in terms of many-path integrals 
with counting intersections (raprochements) of paths (like it was considered in [14] with 
respect to single copy). Under the weak short-correlated bath noise approximation, evolution 
equations are time-local and then, at least, right-hand side of equation for (^(t))^ can be 
represented in standard Lindblad form for generators of quantum contracting semigroups 
(although in our case this is not natural form). But the Langevin equation (20), which 
in our case realizes stochastic dilation of a dissipative semigroup behaviour, lies rather far 
from Langevin equations under use in modern "quantum stochastic calculus" [13] or in well 
developed bosonic stochastic calculus [15]. The peculiarity of Eq.20 is that it takes into 
account phase volume exchange between dynamical subsystem (electrons) and thermal bath 
and, therefore, the dilation it gives is not unitary in momentary sense although eventually 
(after complete averaging) works as unitary one. To the best of our knowledge, such sort of 
Langevin equations was still unknown. If it is so then Eq.20 itself also can be qualified as 
result of this paper. 

In the frame of this equation, we came to very intriguing task about spectral and other 
properties of the Fokker-Planck operator Aq = Afree+ A^.^.+ A^y whose peculiarities are non- 
self-adjointness, non-ellipticity (induced by that of A^^, ) and continuous infinitness of Ker 
Aq . Particularly, Fokker-Planck approximation turns into exact theory if throw out thermal 
bath at all but instead include delta-correlated x'(t) - y'{t) noise equivalent to pair electron- 
electron interactions. In this case Axx = and "bad" operator A^y becomes responsible for 
everything (perhaply, in particular, for such "bad" phenomenon as flicker noise). 

I acknowledge Dr. I.Krasnyuk and Dr. Yu.Medvedev for useful conversations. 
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